A VIABLE COSMOLOGY WITH A SCALAR FIELD COUPLED TO THE TRACE 

OF THE STRESS-TENSOR 
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We study the cosmological evolution of a scalar field that couples to the trace T — T% of energy 
momentum tensor of all the fields (including itself). In the case of a shallow exponential potential, 
the presence of coupling to the trace T in the field equation makes the energy density of the scalar 
field decrease faster thereby hastening the commencement of radiation domination. This effect 
gradually diminishes at later epochs allowing the scalar field to dominate the energy density again. 
We interpret this phase as the current epoch of cosmic acceleration with 57^ = 0.7. A variant of this 
model can lead to accelerated expansion at the present epoch followed by a a(t) oc i 2//3 behaviour 
as t — » 00, making the model free from future event horizon. The main features of the model 
are independent of initial conditions. However, fine tuning of parameters is necessary for viable 
evolution. 

PACS numbers: 



I. INTRODUCTION 

The interpretation of current observations in terms of 
an accelerating universe requires an exotic form of energy 
density with [p + 3p) < 0. A wide variety of scalar field 
models have been conjectured for this purpose including 
quintessence pj, K-essence0, tachyonic scalar fields 0,0, 
[3 with the last one being originally motivated by string 
theoretical ideas 0. (For a recent review of issues related 
to the cosmological constant and dark energy, see 0). 
Since all these models have a free, undetermined function 
V(<j)), it is possible to incorporate any reasonable a(t) in 
any of these models 4]. Thus, at a fundamental level, 
all these models can be objected to — quite correctly 
— as lacking in predictive power. The appropriate way 
of judging these models, therefore, will be based on the 
physical input that is used in constructing the model and 
the generic features which arise in the model. 

We investigate here a model in which the coupling of 
the scalar field to the matter is based on a physical prin- 
ciple which can be motivated as follows: The most ob- 
vious generalization of the Newtonian gravity, based on 
V 2 4> = A-rtGp will be U<j> = 4irGT in which the trace 
T = T" of the energy momentum of fields act as a source 
for a scalar field. This theory, however, has to be non- 
linear, since the scalar field itself has a nonzero trace 
for the energy momentum tensor which should appear 
on the right hand side. Further, since cosmological con- 
stant has a nonzero trace for stress tensor, this scalar 
field will also couple to the cosmological constant and 
the effective cosmological constant will become dynami- 
cal. Such a model was developed and investigated more 
than a decade back 0, , as a possible solution to cos- 



mological constant problem. This work, however, was 
ambitious in the sense that it attempted to tackle the 
problem without any fine tuning and did not succeed 0] • 
The current thinking in cosmology is remarkably toler- 
ant to fine tuning of parameters and virtually every scalar 
field model suggested in the literature has either explicit 
or implicit fine tuning of parameters. In view of this, one 
can resurrect the above model, in which the form of the 
interaction arises from a clear physical requirement, and 
explore its consequences for cosmology. We attempt to 
do this in this paper and show the the resulting models 
have very nice properties. 



II. DYNAMICS OF THE SCALAR FIELD WITH 
A COUPLING TO T a a 

The Lagrangian for a system made of the gravitational 
field, matter sources and the scalar field which couples to 
the trace of the total energy momentum tensor is derived 
in [j| . (This derivation is briefly summarized in Appendix 
A for completeness.) The complete action, which takes 
into account the trace coupling of the scalar field, has the 
form (see Appendix A): 



S = {WttG)- 1 J Ry/=g~d i x+^ J f3{(j))(t) l (j) i 



(8ttG)- 



*{4>)V{cj>)yf=g-d A x + S s , 



-gd x 



(1) 



where the coupling to the trace is given by an arbitrary 
function f(<j>) and 



a(cf>) 



1 



1 + 4/(0)' 



1 



1 + 2/(0) 



(2) 
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The S'sourco is the remaining source term (radiation, mat- 
ter ... ), with appropriate coupling to (see Appendix 
A). The energy momentum tensor for the field which 
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arises from the action is given by 

r k = a^)v{<t>)g ik + m 



V4> k - \g lk 4>>4>3 



(3) 



We shall assume that scalar field is evolving in an 
isotropic and homogeneous space time and that <j> is a 
function of time alone. The energy density p^ and pres- 
sure pc/, obtained from T tk are 



P4> = 0{<t>) y + a(4>)V(4>), P4> = p{ft>)1j - a(<l>)V(<t>) 



The coupling to T" modifies the field energy density 
pressure p^ and the equation of state parameter 



(4) 



In a universe containing dust like matter and radiation, 
the action (JIJ leads to following evolution equations for 
the field: 



3a 
— < 



1 + 2/ 



V+4/f ' 



v ^ ( 1 + 2 /) _i_ (L f (i + v) 



(5) 



(The radiation, being traceless, does not couple to <j>.) 
The Friedman equation with the modified energy mo- 
mentum tensor given by J2J) is 



a 2 (t) _ 8ttG 



<P V{4>) 

2(1 + 2/) 1 + 4/ 



Pb{a) 



(6) 



where background energy density due to radiation and 
matter is given by 



with 



C(«) 



1 



(K e /P e ) + 1 



where the ratio of effective kinetic to potential 
energy(X e /P e ) is given by 



Ke 

p. 



«(</>) 2V(cf>) 



(11) 



Interaction modifies the kinetic as well as the potential 
energy through a(4>) and f3(4>), however, p^ bears the 
same relation to their ratio as in the case of standard 
scalar field cosmology. Equation ©, of course, is not 
valid in presence of matter; neverthless numerical work 
shows that the ratio (K e /P e ) is still a good parameter 
which influences the dynamics. The evolution of poten- 
tial to kinetic energy ratio plays a significant role in the 
growth or decay of the energy density at a given epoch 
and will be crucial in the following discussion. 

Let us next address the question of choice of potential 
of the scalar field which would lead to a viable cosmolog- 
ical model. The obvious restriction on the evolution is 
that, starting from Planck's time, the scalar field should 
survive till today (to account for the observed late time 
accelerated expansion) without interfering with the nu- 
cleosynthesis of the standard model. 

Exponential potentials lead to solutions which are the 
attractors of evolution equations when g c = and pro- 
vide backdrop for the understanding of dynamics in our 
case. A standard model with 



Pb(a) 



PR , Pm 

a 4 + a 3 (l + /) 



(7) 



The simplest form of coupling which we shall adopt is 
the linear one 



/(</>) = 9 cWM p ) 



(8) 



where M p = y/l/(8ivG) is the reduced Planck's mass and 
g c is a coupling constant. Since equations (J5J and © re- 
duces to the standard form for g c = 0, this parameter 
gives the strength of the forcing term. Equation © de- 
scribes the evolution of the scalar field in the expanding 
universe under the influence of an external "force" , which 
depends upon the field and its kinetic energy (l/2)0 2 . 
The influence of this term on the dynamics of the scalar 
field is accumulative. 

In the absence of matter, the conservation equation 
formally equivalent to 0, has the usual form 

p<f, + 3H(p<f > +p < f,) = (9) 

and the evolution of energy density is given by 

- f 6(l-C(o))— rm\ 

Pd> = A^e J ^ " a (10) 



V(<t>) = V exp 



(12) 



(and g c = 0) is well studied in literature and can be di- 
vided into two categories: (a) The exponential potentials 
for which p^ scales slower than the background density 
Pb = (\/a n ) which translates to A < y/n. (b) Potentials 
for which scalar field energy density scales faster than 
the i.e. A > y/n. In the first case if the background 
energy density was sub-dominant, it would become more 
sub-dominant and radiation domination will never occur. 
In the second case there exists a scaling solution which 
mimics the background energy density that is dominant, 
with 



p<t> 



P4> + Pb 



n 

a7 ; 



i 

P0 °c — 

a n 



(13) 



The discussion of this model with g c ^ (so that the 
coupling to the trace of the stress tensor is switched on) 
is given in Appendix B. It turns out that this model is 
not as attractive as another variant of the potential with 

V(4>) = Voe M p which we shall concentrate on next. 



3 



III. QUINTESSENCE WITH EXPONENTIAL 
POTENTIAL IN THE PRESENCE OF THE 
COUPLING TO THE TRACE T 

We shall now consider a model with 

V(<j>) = V exp(A^), (A < V2) (14) 

in which the field evolves towards the origin rather than 
away from the origin. In order to investigate the dynam- 
ics described by (J5J and ©, it would be convenient to 
cast these equations as a system of first order equations 



Y{ = Y 2 



(15) 



Y' 

1 2 



-Y4 



9c 



H(Y U Y 2 ) L(l + 2 5c y x )' 2 
(1 + 2g c Y 1 ) - dViYx) (1 + 2g c Y 1 ) 



(l + 4 5c Yi 
fa (1 + 2g c Yy 
a 3 



V{Yi) 



dYx (l + 4ff c y x ) 



(i + g c Y 1 y 



(16) 



where 



Yi = — y? = — — 



V 



V(Yi. 



(17) 



and prime denotes the derivative with respect to the 
able N = ln(et). The function H(Y\,Y2) is given: 



\ 



Y 2 



V(Y 1 



Pb 



M4 



2(l + 2 9c F0 (l + 4 5c Yi) 

(18) 

where p b = p l r e~ m + p l m e- 3N /{l + g c Yi). The initial 
conditions will be set at Planck time with a% = 1. (The 
scale factor today would nearly be 10 31 .) Further, for 
the background energy density, p r = p l r /a , we make the 
choice, p\ — based on the point of view that, at the 
radiation dominated phase near the Planck's epoch, the 
energy density of radiation was of the order of Planck en- 
ergy as it was the only scale available at that time. The 
initial values for (j> and <\> as well as the the parameters 
in the potential are chosen so as to ensure a viable cos- 
mological model. We take the initial values Y\ = 550; 
Y-2 = and A = 0.52, g c — 5.9. The value of Vq is chosen 
so that initially, p\/ ' p % R ~ 10 2 , Of these choices, A and g c 
are important for obtaining a viable model; Y\ is arbi- 
trary except for the condition Y\ 1. The choices made 
for Yi and Vq are merely done for the sake of convenience; 
different choices will lead to same qualitative behaviour. 
Some of these choices could be construed as fine tuning 
in the model. But it is no worse than the fine tuning in 
any other scalar field model. 

The model described by (|14|l supports an eternal in- 
flation in absence of the coupling to the trace T . The 
forcing term induced by the back reaction in the field 



evolution equation acts against the Hubble damping and 
hastens the evolution of field energy density p^. We have 
plotted the ratio of kinetic to potential energy of the field 
in figure ^ It is seen that due to coupling to the trace 
T , the kinetic energy acquires a large value making the 
potential term unimportant and pushing the field into ki- 
netic regime. As a result, the field energy density crosses 
the background energy density and continues evolving 
as (1/a 6 ) for quite some time (see figure [SJ. After the 
crossover, the dominant contribution to the Friedman 
equation Q comes from the background energy density 
Pb- This in turn leads to the increase in damping and con- 
sequently slowing down of the decay of p^ allowing the 
ratio / Pb) to grow. The energy density gets locked to 
a constant value in this process and acts like cosmological 
constant thereby making to grow. The forcing term 
containing 4> 2 in the field evolution equation is ineffective 
in this regime. 

A remarkable thing happens when the field starts 
rolling again at the end of locking regime. As discussed 
above, the ratio of kinetic to potential energy of the field 
(K e /P e ) provides a yardstick to determine the energy 
scaling. At the end of locking period, the field again 
starts rolling down the hill. Though back reaction again 
builds up, the kinetic to potential energy ratio is now 
much smaller than it was in the beginning when was 
dominant. The ratio first increases, then fluctuates near 
its background value, finally heading towards zero (figure 
EJ. As a result, p^ first (nearly) tracks the pb, gradually 
approaches it and ultimately overtakes it. This allows 
the scalar field to dominate and gives rise to accelerated 
expansion at late times. By suitably tuning the param- 
eters of the model, it is possible to get fi^ = 0.7 today 
(figure 0J . In figure we have displayed the equation 
of state parameter w<j, to demonstrate the convergence 
of different initial conditions to the inflationary attrac- 
tor. Figure © shows that the w$ mimics de-Sitter like 
behavior shortly after the current epoch. 

The role of tp 2 term on the right-hand-side of evolution 
equation (j^J is very important in the beginning. The 
size of the undershoot crucially depends upon the cou- 
pling constant g s . Larger value of g s leads to deeper 
undershoot because the ratio K e /P e is sensitive to cou- 
pling initially. After the locking period this ratio of the 
order of few as background (radiation) dominates over 
the field this time. During evolution, the scalar field 
moves from larger values towards the origin. As a re- 
sult the influence of the back reaction in the Friedman 
equation © decreases and the field energy density grad- 
ually moves towards the background, overtakes it and 
joins with the (non-coupled) inflationary attractor(This 
may be contrasted with the behaviour in the model de- 
scribed by H12JI . in which the scalar field evolves from 
origin to larger values. In this case the contribution of 
field energy density in the Friedman equation gets gradu- 
ally suppressed allowing the background to dominate for 
ever. In this model tracks the background and never 
comes to dominate it(see figures. H3I14I and ll5H . This de- 
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velopment is displayed in figure (JTJ which clearly shows: 
(i) deeper undershoot for larger coupling and (ii) joining 
of the inflationary track earlier for smaller value of g s . 
Figure [7| also shows that the field continues to stay in 
the usual (non-coupled) attractor phase for quite some 
time till the potential terms induced by the coupling in 
field evolution equation become important. This happens 
at very late times and the influence of these terms will 
be discussed below. As for the matter term, numerical 
analysis shows that they do not influence the dynam- 
ics significantly. Thus coupling plays an important role 
in the intermediate regime (which also maifests in figures 
12151 and ISjl thereby allowing to produce current observed 
accelerated expansion with necessary tuning of g c along 
with a. 

Let us briefly summarize the role of the coupling to the 
trace T .The exponential potential under consideration is 
shallow and supports an eternal inflation in the standard 
scenario. The introduction of the coupling to the trace T 
hastens the decay of p<f>, pushing the scalar field into ki- 
netic regime. After the locking period, its role gradually 
diminishes making it possible for the scalar field, at late 
times, to join back the inflationary track it was detracted 
from. At this stage, the scalar field is approaching the 
origin and we would expect the equation of state w<p to 
settle in its attractor value. However, the equation of 
state becomes equal to (—1) and stays there — which is 
evident from figure The same effect is reflected in the 
behavior of p^ shown in figured where the field energy 
density becomes constant after staying in its attractor 
value for some time. This seems to be a general feature 
of the model under consideration and requires explana- 
tion. 

The coupling to the trace T modifies the field evolution 
through three extra terms on the right hand side of (16). 
The first term contains cb 2 and plays important at the 
early epochs and during tracking as discussed above. The 
second term contains the potential which contributes in 
the potential dominated regimes: during the turn around 
of and when the field is rolling slowly. We will show 
that this term is responsible for the cosmological constant 
like behavior of the scalar field at very late times , i.e., 
a — > oo. In order to see this, we will cast (|15H in the form 



% + a -Y 2 
a 



Y- 



-VtotiYi 



1 + 2g c Y 1 dYi 
where Vtot(ii) is the effective potential given by 



(19) 



Hot(ii: 



(1 + 2g c Y 1 ) 



(l + 2.g c ri) 



V(Y 1 )dY 1 



(l + 4 flc yi) 3t (l + 4ff c Fi) 2 

(20) 

where V = (Vo/M^e 1 . (The matter term is sub- 
dominant in this limit and is hence ignored.) In the 
model described by (JT^J, the scalar field is evolving from 
large values towards the origin and can take negative 
values. But when Y\ — » —(l/2g c ) the system will be- 
come singular because of the divergence of the first term 
in i|19|) . However, the second term in the expression of 



Vtot(Y\) does not let it happen. In fact, this term sets 
an infinite potential barrier at Y\ = — (l/4g c ) and does 
not allow the field to reach the singularity. The potential 
which was attractive away from the origin becomes now 
repulsive followed by the infinite barrier on its left (see 
figure IHJ). This happens when 



Y c 



4:g c - A 
4Xg c 



(21) 



At this stage, the field <j> is rolling slowly. It crosses the 
critical point and moves towards the barrier but due to 
insufficient kinetic energy gets reflected back by the bar- 
rier. After making few such oscillations, the field settles 
at the critical point permanently and mimics the cosmo- 
logical constant. For the model described by (|14fl with 
A = 0.52 and g c = 5.9, Yf = (j) cr /M p ~ 1.88 which is 
confirmed by numerical integration (see, figure 



IV. MODELS OF QUINTESSENCE FREE 
FROM FUTURE EVENT HORIZON PROBLEM 

In the model described above, like in any generic model 
of quintessence, universe will be accelerating forever lead- 
ing to future event horizon — which could pose a problem 
in the string theoretical context. It is, however, possible 
to tackle this problem by using any potential which is 
exponential for large <f> and is a power law near the ori- 
gin: That is, we require V(<fi) — > exp(A0/Mp) for large 
(j) and V{<t>) -> {Xcf/MpfP for -> 0. The power law 
behaviour near the origin will lead to oscillations when 
the field approaches <j> ~ 0. The mean equation of state 
of the scalar field is then given by ^3 



V{cj>)/ 



P-1 
p + 1 



(22) 



As a result the scalar field energy density and scale factor 
have the following behavior 



-3(l + <u») 



ffl0ct |(i+<«»»- 



The scalar field behaves like pressure less dust for p = 1 . 
By suitably adjusting the parameters in the model one 
could ensure that the oscillations occur in the future, i.e., 
well after the present epoch when f2^ has reached value 
equal to 0.7. One can avoid the future event horizon in 
such a model, irrespective of the explicit form of V{(j>) 
which is used, as long as it interpolates between an expo- 
nential and power law behaviour. One of the many pos- 
sible choices of interpolation is a potential of the form 
V{4>) = Vbsinh 2p (A</>/2Mp), for which we have evolved 
the field equations numerically. The results for a partic- 
ular choice of parameters are shown in figures l|10|) and 

We summarize below the chronology of main events in 
any such model with a shallow potential which behaves 
like a power law near the origin: (i) The scalar field starts 
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rolling slowly in the shallow potential. The coupling to 
the trace T almost instantaneously kills any inflationary 
behavior and makes the scalar field to enter into the ki- 
netic regime. The scalar field continues to stay in the 
kinetic regime (with p^ oc a~ 6 ) and soon p^ drops below 
the pb- (ii) Once this happens, the dominant contribution 
to the Friedmann equation comes from the background. 
This, in turn, allows p^ to turn towards pb leading to a 
locking regime during which the value of p^ stabilizes and 
remains relatively unchanged for a considerable length of 
time and the scalar field behaves like cosmological con- 
stant, (iii) The field starts rolling down the potential 
again and the back reaction builds up such that the ki- 
netic energy of the field exceeds its potential energy by 
a factor of few in the beginning. However, as the field 
rolls towards smaller and smaller values, the role of back 
reaction gradually diminishes and the potential energy 
starts catching up again with the kinetic energy. As a re- 
sult p^ slowly approaches the back ground energy density 
and finally overtakes it. The scalar field is now back on 
the original inflationary track. This phase is interpreted 
as the late time accelerated expansion of the universe, 
leading to present day value of = 0.7. (iv) After the 
current epoch, the universe continues in the phase of cos- 
mic acceleration for quite some time till the scalar field 
approaches <p ~ and starts oscillating about it giving 
rise to < w,p >= and turning the scalar field to behave 
like cold dark matter (see figure ITT)> . This is expected 
to happen in future and the model is free from the event 
horizon problem at late times. (There are other possible 
solutions to this problem; see, for example, @, El). 



V. CONCLUSIONS 

We have examined the dynamics of a scalar field which 
couples to the trace of all the field including itself in FRW 
background. The coupling to the trace T of the field mod- 
ifies the energy momentum tensor and induces a forcing 
term in the field evolution equation. The forcing term 
acts against the Hubble damping changing the evolution 
of field energy density significantly. In case of a steep 
potential, it amounts to increasing the steepness of the 
potential. In case of the shallow potentials in which the 
field rolls slowly and mimics cosmological constant, the 
coupling to the trace T plays an important role. It forces 
the scalar field into kinetic regime allowing the radiation 
domination to commence. We have discussed a model in 
which the effect of the back reaction gradually diminishes 
allowing the scalar field to catch up with the inflation- 
ary attractor regime at late times, thereby accounting for 
the observed cosmic acceleration with 17 ^ = 0.7. We have 
also presented a variant of the model in which the scalar 
field plays the role of quintessence at the current epoch 
and mimics dust like dark matter in future which elim- 
inates the future event horizon. The general features of 
the model are shown to be independent of the initial con- 
ditions. However, fine tuning of some of the parameters 



is necessary to achieve viable evolution. 

APPENDIX A: ACTION FOR SCALAR FIELD 
COUPLED TO THE TRACE OF STRESS 
TENSOR 

We summarize here the derivation of the self-consistent 
Lagrangian for a scalar field coupled to the trace of the 
energy momentum tensor of all matter, originally given 
in |3, for the sake of completeness. 

Consider a system consisting of the gravitational field 
g a b, radiation fields, and a scalar field <f> which couples 
to the trace of the energy-momentum tensor of all fields, 
including its own. The zeroth order action for this system 
is given by 

(Al) 

where 

A grav = {WnG)- 1 [ R^f~gatx- [ A^d 4 x, (A2) 



A f = \j WiV^gd^x; A ( n°t = v I TfWMV^g** 

(A3) 

Here, we have explicitly included the cosmological con- 
stant term and r\ is a dimensionless number which 
'switches on' the interaction. In the zeroth order action, 
T represents the trace of all fields other than tp. Since 
the radiation field is traceless, the only zeroth-order con- 
tribution to T comes from the A term, so that we have 
T = 4A. The coupling to the trace is through a function 
/ of the scalar field, and one can consider various possi- 
bilities for this function. The constant </>o converts <f> to a 
dimensionless variable, and is introduced for dimensional 
convenience. 

Since the stress-tensor of the scalar field has a non-zero 
trace, we must add to T the contribution T$ — —<$4>i of 
the scalar field. If we now add to T in the interaction 

term A^ t further modifies T^ k . This again changes T$. 
Thus to arrive at the correct action an infinite iteration 
will have to be performed and the complete action can 
be obtained by summing up all the terms (see 9J). The 
full action can be found, more simply, by a consistency 
argument. 

Since the effect of the iteration is to modify the expres- 
sion for A^ and Aa, we consider the following ansatz for 
the full action: 

A = J RV~gd A x- J a(<f>)A^d 4 x 

+ \ J (3(cbW&V=gd 4 x + A Ilid (A4) 

Here oc((j>) and /?(</>) are functions of </> to be determined 
by the consistency requirement that they represent the 
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result of the iteration of the interaction term. (Since 
radiation makes no contribution to T, we expect A rac j to 
remain unchanged.) The energy-momentum tensor for <p 
and A is now given by 



so that the total trace is T to t = 4a(^)A — /3(0)</> ! 0;. The 
functions ot{4>) and /?(</>) can now be determined by the 
consistency requirement 



V 



(A5) 



- I Ay^gd x + - I 4? 4>i\/^g d x + rj \ T tot f {cj> / O ) y/^g d x 



(A6) 



r 



Using T to t and comparing terms in the above equation 
we find that 



a(4>) = [l + 4r)f}-\ /?(</») = [1 + 277/]- 
Thus the complete action can be written as 

i r a 



.4 



16ttG 
1 

+ 2 



Ry/—g d x 



1 + 417/ 



I + 277/ 



f=g d*x + A T 



ad 



(A7) 



-g d A x 



(A8) 



(The same action would have been obtained if one uses 
the iteration procedure.) It is obvious that the method 
works for any source for which the trace of stress tensor, 



T, is proportional to the Lagrangian L. If T = pL, then 
the coupling replaces L by L(l — p,rjf)~ l . [For the kinetic 
energy term = (l/2)(/> a </> Q , we have, T$ = -(f) a <j> a = 
(—2)L ! p and for the cosmological constant term Ta = 
4A = (-4)L A . These lead to the factor (1 + 2-qf)- 1 



and (1 + 477/)" 



.] In general, if the action is a 



homogeneous function of degree D in g ab , the T = (2D)L 
and the coupling factor is (1 — 2Drjf)~ 1 . For a system of 
dust like particles with the action for each particle being 



the integral of ds = \J g a bdx a dx b , we have D = 
and the coupling factor becomes (1 + 77/) _1 (see 



9] 



1/2) 
)• 



In the presence of cosmological constant and sources 
which are conformally invariant, the action in (| A8|) leads 
to the following field equations, 



J 



R 



ik 



1 



-8ttG 



ik ik 



A 
8^G 



T, 



tracclcss 



(A9) 



r 



2 m ' 



A q'Q) 
8^G P((p) 



= 



(A10) 



Here, □ stands for a covariant d'Lambertian, x^ acclcs = 1S 
the stress tensor of all fields with traceless stress tensor 
and a prime denotes differentiation with respect to <p. 
In the cosmological context, this reduces to 



a 1 + 277/ 



A f(l + 2ry/) 
V 2nG (I + 477/) 2 [AU) 



8ttG 



1 



2 1 



A 



1 



Pb{a) 



277/ 87rG(l + 47y/) 

(A12) 

where Pb(a) = p l R /a A . It is obvious that the effective 
cosmological constant can decrease if / increases in an 
expanding universe. The result can be easily generalized 



for a scalar field with a potential by replacing A by V((f>) 
in the action. Similarly, the presence of dust like matter 
adds the term (p^/a 3 )(l + rjf)^ 1 to the source term in 
(|A12(1 and a corresponding derivative term in (|A11|I . 



APPENDIX B: COSMOLOGY WITH 

V oc expt-A^/Mp) AND COUPLING TO T 

In this case, the field evolves to large <j) a t late times 
making the matter term a~ 3 [l + 77/] _1 sub-dominant to 
the radiation term a~ 4 . We will, therefore, ignore the 
matter term for simplicity. We first discuss the possi- 
bility when A > y/n. In this case the potential is steep 
leading to p^ cx 1/a 6 . This makes p^ sub-dominant rela- 
tive to pi, and the motion will become strongly damped 
allowing the kinetic energy to decrease. This, in turn, al- 



7 



lows to catch up with the background and stay there 
for ever till the attr actor regime is reached. The role of 
coupling to T will be to increase the kinetic energy and 
keep p0 in the kinetic regime for longer. The will catch 
up again with the background, though it takes longer to 
reach the attractor (scaling) solution with the duration 
depends upon g c . Therefore, in this case, the role of the 
forcing term is equivalent to increasing the steepness of 
the potential given by A (see figure lL3"j) . Thus nothing 
new or interesting happens in this case . 

The case with A < ^Jn ( n — 4 in our case) is more 
interesting. We shall be working with A < \[2 in which 
case the exponential potential supports a never ending 
inflation. The forcing term which acts against damp- 
ing. Due to coupling to the trace T the kinetic energy 
acquires a large value making the potential term unim- 
portant and pushing the field into kinetic regime. De- 
pending upon the value of coupling g c it may take longer 
or shorter to force p§ to come into kinetic regime which 
then crosses the background energy density, (see figure 
I15fl . The extent of the overshoot depends upon the cou- 
pling constant g c . The evolution of p$ in this model pro- 



ceeds more ore less in the similar fashion as in case of the 
model described by ljl"4"|) till the turn around. At the end 
of locking period, the field again starts rolling down the 
potential. Though back reaction again builds up, the ki- 
netic to potential energy ratio is now much smaller than 
it was in the beginning when p^ was dominant ffigurc lT4l . 
The field does not join the inflationary attractor in the 
present case. Because as it runs down the hill away from 
the origin to larger and larger values, the contribution 
of V(4>) and 4> 2 /2 to the right hand side of Fricdmann 
equation © gets suppressed systematically and p^ re- 
mains sub-dominant to the background (radiation). The 
ratio (K e /P e ) settles at value of 2 ( = 1/3) and the 
scalar field tracks the background and gets locked there 
for ever. 
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FIG. 1: The ratio of kinetic to potential energy of scalar 
field is depicted for the exponential potential 1141 ) with A = 
5.2 , 1/ 1/4 ~ 2.8 x l(T 30 Mp with two different values for 
coupling (chosen for illustration): g c — 0.04 (dashed line) 
and g c = 0.041 (solid line). The kinetic energy ( as compared 
to the potential energy ) is seen to build up fast to a large 
value forcing p^ into the kinetic regime which otherwise would 
scale slowly with the scale factor (p^ oc a ) for the model 
with the potential in 1141 . 
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FIG. 2: The energy density is plotted against the scale factor: 
solid line corresponds to p^ for g c — 5.9 in case of the model 
described by JHJ with A = 0.52 and V ( 1/4 ~ 2.8 x 10" 30 M P . 
The dashed and dotted lines corresponds to energy densities 
of matter and radiation respectively. The large value of the 
scalar field kinetic energy (relative to the potential energy ) 
induced by the forcing term in the evolution equation ensures 
that the p^ drops below pt, after which remains approxi- 
mately constant for a period of time. At the end of this phase, 
the scalar field tracks the background energy density before 
becoming dominant and driving the current accelerated ex- 
pansion of the universe. 
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FIG. 3: Evolution of the kinetic to potential energy ratio of 
scalar field for the model described in figure JSJ starting from 
the locking regime onwards. At the end of locking period, the 
ratio fluctuates about its background value and then nearly 
tracks it slowly decreasing towards its attractor value. The 
ratio does not lock itself to its attractor value and approaches 
zero where it stays asymptotically. The vertical dashed line 
marks the current epoch. 
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FIG. 4: Evolution of dimensionless density parameter Q is 
shown as a function of the scale factor for the model described 
in figure for: (i) the scalar field (solid line), (ii) radiation 
(dotted line) and matter (dashed line). Late time behavior of 
the scalar field leads to the present day value of 0,$ = 0.7 and 
On = 0.3. 
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FIG. 5: The equation of state parameter of scalar field for the 
exponential potential 11411 showing different initial conditions 
converge to attractor solution at late times. The horizontal 
dashed line depicts the = 1/3 = wr. The solid and dashed 
lines start with widely different initial conditions but lead to 
similar behaviour at late times. 




FIG. 6: Evolution of equation of state parameter of scalar 
field. The coupling to the trace T brings scalar field (which 
was slowly rolling down the shallow potential) into kinetic 
regime. The scalar field stays in this regime for a long time 
before p^ grows towards pb- During this phase, w$ gets locked 
to lo^ = — 1 ( locking regime). After the locking period ends, 
pc/y tracks the background energy density slowly approaching 
it from below with equation of state w$ ~ 1/3. At late times 
overtakes the background and becomes dominant to ac- 
count for the current epoch of cosmic acceleration. During 
this time decreases and heads towards its inflationary at- 
tractor value (wtf, = (A 2 — 3)/3)). However, before it could 
settle there, the slowly moving field <fi towards the origin gets 
trapped near the barrier formed there due to coupling to the 
trace T . The equation of state locks to —1 permanently. The 
vertical dashed line marks the present epoch. 
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FIG. 7: The evolution of energy density in the field is shown 
for two values of g c — 0.9(solid line) and g c = 0.07(dotted 
line) for the model described ll 1 1 li with A = 0.52. The ra- 
diation energy density is depicted by the dashed line. The 
second dashed line below shows the(non-coupled) inflation- 
ary behavior of p^>(p^ oc a~ x ). The decay of p<^ slows down 
earlier for smaller value of the coupling constant g c and ap- 
proaches the inflationary attractor behavior after overtaking 
the background (radiation) energy density. After staying with 
the attractor value for quite some time, p,/, becomes constant. 
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FIG. 8: The effective potential in equation (120ft for A = 
0.52, g c = 5.9. This effective potential becomes repulsive near 
the origin when (<fi C r/Mp) = 1.88 followed by an infinite bar- 
rier at (<j> /M p ) — —0.042 which shields the singularity situ- 
ated at ((j) s /M p ) = -0.084. 
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FIG. 9: Late time evolution of the scalar field (solid line). Af- 
ter the current epoch, the scalar field passes across the critical 
point, moves towards the barrier and gets reflected in the ef- 
fective potential shown in the previous figure. It makes few 
oscillations and settles at the critical point (<f) cr /M v ) = 1.88. 
The dashed line depicts the inflationary attractor behaviour of 
the scalar field in case of the exponential potential 1141 with 
A = 0.52 (4>/M p — — Aln(a/ai) + constant) and the dotted 
line corresponds to {4> cr /Mp) = 1.88. 
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FIG. 10: The dimensionless density parameter is plotted as 
a function of the scale factor for: (i) scalar field (solid line) 
with a potential V = Vo s'mb 2p (\<f>/2M p ), radiation ( dotted 
line ) and matter (dashed line) . The late time behavior of the 
scalar field is shown to lead to the present values of Q$ — 0.7 
and Q, m — 0.3 
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FIG. 11: Evolution of energy density is shown for scalar field 
(solid line) radiation (dotted) and matter (dashed) for a po- 
tential V oc sinh 2p (A<^/2M p ). The energy density drops 
below pbg, turns around catching up with it in the matter 
dominated era. The scalar field dominates at late times and 
derives the current cosmic acceleration. The universe contin- 
ues in this phase for some time till field oscillations ( near 
the origin <j> = ) build up in the system leading to average 
equation of state < >= 0. The field behaves as cold dark 
matter thereafter. 
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FIG. 12: Evolution of the scalar field is shown for the model 
described in figure [TT1 in the asymptotic future. (The current 
epoch is a ~ 10 31 .) After driving the accelerated expansion 
for quite some time, the field enters into oscillatory phase near 
the origin with the average equation of state < >= 0. 
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FIG. 13: The evolution of scalar field density corresponding 
to g c — 0.05 (solid line) and g c — (dotted line) is shown for 
the potential 1121 with A = 6. The radiation energy density is 
shown by the dashed line. The pattern of tracking shows that 
the effect of forcing terms amounts to increase the steepness 
of the potential. 
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FIG. 14: Evolution of the kinetic to potential energy ra- 
tio of scalar field for the potential (1 1 2t with A = 1.2 and 
g c — 1 starting from the locking regime onwards. At the end 
of locking period, the ratio fluctuates about its background 
(radiation in this case) value (K e /P e = 2) and ultimately 
approaches it. 
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FIG. 15: The energy density of scalar field p$ is plotted for 
g c = 1 (solid line) and g c =0.1 (dotted line) for shallow ex- 
ponential potential 1121 with A = 1.2. The dashed line shows 
radiation energy density. The effect of the forcing term in the 
evolution equations is shown to accelerate the scaling of p^,. 
The stronger coupling is seen to hasten the commencement 
of kinetic regime. The scalar field ultimately tracks the back 
ground and gets locked up there for ever. 



